Resummation of jet mass with and without a jet veto 
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Calculating the distribution of jet masses in high-energy collisions is challenging because fixed- 
order perturbation theory breaks down near the peak region, and because multiple scales complicate 
the resummation. To avoid using a jet veto, one can consider inclusive observables, in which every 
particle is in a jet. We demonstrate that calculating the mass of the hardest jet in multijet events 
can be problematic, and we give an example of an inclusive observable, asymmetric thrust, which 
can be resummed to next-to-next-to-leading logarithmic accuracy. Exclusive observables with out- 
of-jet regions are more complicated. Even for e + e~ dijet events at energy Q, to calculate the mass 
m of jets of size R, one must impose a veto on the energy u of extra jets to force dijet kinematics; 
then there are both logm/Q and log m/ui singularities. To proceed, we suggest a refactorization of 
the soft function in the small R limit To justify this refactorization, we show that the expansion of 
the resummed distribution is in excellent agreement with fixed order. This motivates considering 
the expansion around small R as a useful handle on producing phenomenologically useful resummed 
jet mass distributions. The strong evidence we give for refactorization at small R is independent of 
non-global logarithms, which are not the subject of this paper. 
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1. INTRODUCTION 

The Large Hadron Collider has already produced enor- 
mous numbers of high energy jets. These jets provide 
a wealth of information about QCD as well as a vital 
area for exploration of new physics. While much activity 
has been devoted to the calculation of the distribution 
of jets, for example using leading order (LO) or next-to- 
leading order (NLO) perturbative calculations in QCD, 
very little has been said with any precision about the 
substructure of the jets themselves. This is unfortunate, 
because the substructure of jets may be critical to find- 
ing new physics. In fact, there has been much progress 
over the last few years on practical methods using jet 
substructure to separate signal from background at the 
LHC [lri. However, even for the simplest of jet shapes, 
the jet mass, these studies are forced to rely on the Monte 
Carlo approximation. Needless to say, it would be great 
to have a systematically improvable way to calculate jet 
shapes, and it may even prove essential if new physics 
shows up in one of these channels. 

The main theoretical difficulty with calculating sub- 
structure is that fixed-order perturbation theory is in 
general a very bad approximation. For example, at 
LO, jet mass distributions for small mass diverge as 
da ~ 1/ m 2 , while the measured distributions have a peak 
at small mass and then go to zero. The physical peak 
structure is easily explained qualitatively as being due to 
large logarithmic terms a\ log 3 (m 2 /E 2 ), which dominate 
over the fixed-order expansion in a s for small enough 
mass. The leading large logarithms (LL) are resummed 
with the Monte Carlos (or analytically [J]); however, un- 
til recently, it was impossible to go beyond this order. 
The difficulty stems from the fact that while the leading 
Sudakov double logarithms are universal, subleading log- 
arithms are not. The coefficients of the subleading logs 



depend on many variables besides the mass, including the 
jet size, the jet algorithm, and the directions and energies 
of the other jets. Sorting out all these effects and all the 
scales is a daunting task. However, with recent advances 
in effective field theory, it is now possible to make sys- 
tematically improvable calculations of jet substructure, 
at least in some kincmatical regimes. 

While we are mainly interested in jet masses at hadron 
colliders, from a theoretical perspective it is instructive to 
start with e + e~ machines. In e + e~ collisions, jet shapes 
can be approached by first studying hemisphere jets and 
inclusive event shapes like thrust. The factorization the- 
orem in Soft-Collinear Effective Theory (SCET) writes 
the doubly differential hemisphere mass distribution as a 
convolution of a hard function, describing the short dis- 
tance e + e~ — ¥ qq process, jet functions, describing the 
collinear degrees of freedom in the hemisphere jets, and 
a soft function, which describes radiation which can con- 
nect the jets to each other @, Q: 
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Here mr, (run) is the invariant mass of the 4- vector sum 
of the momenta of particles in the left (right) hemisphere, 
with respect to the thrust axis, and Co is the born cross 
section for e + e~ to dijets. The challenge which we ad- 
dress in this paper is to modify this result to describe jets 
of arbitrary size R while still resumming all of the large 
logarithms. 

Before considering non-hemisphere jets, a few com- 
ments are in order about non-global logarithms (NGL). 
Non-global logarithms Q arise in situations where there 
are multiple scales. They become relevant, for example, 
when one integrates the double differential mass distri- 
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bution over m 2 R to produce the left-hemisphere mass dis- them if the distance 
tribution. The full QCD calculation of da/dm 2 L involves 
integrals over hard emissions into the right side which are 
not in the soft or collinear sectors present in the leading 
SCET Lagrangian. To account for these non-global logs 
in SCET would require incorporating 3-jet and higher 
order operators in SCET (see [1,0). These logs appear 
very difficult to resum since there is no known simpli- 
fication for multiple hard emissions in QCD. However, 
they can be avoided by carefully chosen observables. For 
example, thrust and heavy jet mass do not have the non- 
global problems that left-hemisphere or light-jet mass 
have. 

There are other types of non-global logs that are rel- 
evant even for an observable which is well described by 
the leading operator in SCET. For example, consider the 
doubly differential jet mass distribution in the regime 
rriL <C uir <C Q. In this regime, the above factorization 
formula is valid because hemisphere masses (jul and tur) 
are small. This is in contrast to the left-jet mass, where 
in integrating the right-jet mass to its kinematic limit 
one enters the regime with tur ; ~ Q where the factor- 
ization formula fails. So when ra^^mR <C Q, the exact 
distribution is completely dominated by the IR degrees 
of freedom included in SCET. However, SCET does not 
guarantee that the logarithms of and vtlr can be re- 
summed using rcnormalization group evolution (RGE). 
In particular, logs of the form log {rn 2 L /rn R ) are not com- 
pletely fixed by the RGE. If by brute force one could 
calculate to all orders in SCET, the result would repro- 
duce the rriL and mg dependence of full QCD when both 
masses are small 33(. In this paper, we do not attempt 
to predict or resum these challenging NGLs. Instead, we 
give evidence that a certain soft function refactorizes in 
the small R limit. This provides an orthogonal direction 
in the study of jet shapes to the investigation of NGLs. 

Now consider non-hemisphere jets. There has already 
been some work in SCET on understanding finite jet-size 
effects. For example, jet functions for the fcy and cone 
jets have been calculated at NLO [H, [r| , agreement 
with QCD calculations of the 2-jet limit has been shown 
for small R , and RG invariance at the 1-loop level has 
been demonstrated for jet mass and jet angularity distri- 
butions fl3l . [l5j . These works have provided a number of 
non-trivial consistency checks on the SCET formalism. 
However, it has not yet been demonstrated that any of 
the jet shapes exponentiate. Due the presence of multi- 
ple scales, such as the jet size R and the amount of in-jet 
and out-of-jet radiation, there is plenty of opportunity 
for logs like \og(m\/m R ) in the hemisphere case, which 
are not fixed by RGE, to inhibit resummation. 

In this paper, we will consider various jet mass observ- 
ables. We focus on e + e~ — > dijets events, since analyt- 
ical calculations can be performed to order a s and nu- 
merical calculations to a 2 . To lend concreteness to the 
discussion, all events are to be clustered using the Cam- 
bridge/ Aachem (CA) algorithm (see [IH). That is, take 
the pair of particles with the smallest angle and merge 
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is less than R. If no particles have Rij < R then stop. 
The jets can be ordered by energy, and Ei,p^ will denote 
the energy and 4-momentum of the i-th most energetic 
jet. This parameter R can be translated into other jet 
size measures, such as ones used at hadron colliders. 

We begin by considering observables which avoid a jet 
veto. To do this, one has to have an inclusive observ- 
able whose singular region forces the threshold in which 
both jet masses are small. The example we consider is 
asymmetric jet masses and asymmetric thrust, which are 
defined in Section[2l In this case, we show that the asym- 
metric jet mass of the hardest jet does not agree with any 
calculation in SCET, while the asymmetric jet mass of 
the quark (or anti-quark) jet cannot be defined in QCD. 
When one averages over both jets, the results do agree, 
and we check at order a 2 that the next-to-next-to-leading 
logarithms are being resummed. 

Next, we move on to observables which have an out-of 
jet region. To resum the large logarithms of jet mass in 
this case, one must reject events in which there is hard 
radiation in the out-of-jet region. To do this, we intro- 
duce a parameter lu which controls the out-of-jet energy. 
Resumming logs of jet mass in this case is trickier than 
for asymmetric thrust. We argue that expanding around 
a different threshold R ~ 0, provides a new handle on 
these complicated multiscale observables. We postulate 
a refactorized form of the soft function, valid at small R, 
and then show it gives much improved agreement with 
the NLO predictions of QCD, than using SCET alone. 
This indicates than an expansion around small R maybe 
be a promising direction for precision jet substructure 
calculations. 



INCLUSIVE JET MASS: 
THRUST 



ASYMMETRIC 



In this section, we consider a simple inclusive jet mass 
observable, introduced in [13]. Cluster the particles into 
jets of size R using Cambridge/ Aachen (see above). Then 
asymmetric thrust ta is defined as the sum of the mass- 
squared of one primary jet and the mass-squared of ev- 
erything else in the event. Which jet is chosen to be 
primary will be discussed shortly. 

There are three main reasons for considering this odd- 
sounding observable: First, every particle contributes to 
the value of ta- So the region with ta <C 1 forces dijet 
kinematics and the SCET factorization formula should 
hold. Second, there is _R-dcpcndcncc in the jet mass so 
one can study i?-dcpcndcncc of threshold logs. This pro- 
vides a warm up to more complicated exclusive jet mass 
observables. Third, asymmetric thrust may be useful to 
measure jet masses at hadronic colliders. Indeed, it was 
suggested first in the context of dijet events in pp colli- 
sions. But more simply, in direct-photon[l8|] or W/Z [l9[ 
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events at high pt , asymmetric thrust would measure the 
same singular behavior of the jet mass as any other jet 
mass measure, but it may be easier to compute and is 
likely to undergo dynamic threshold enhancement. The 
hadron collider applications are not the subject of this 
paper. 

In dijet events at e + e~ machines, the first question 
is which jet should be the primary jet for asymmetric 
thrust. In data (or in full QCD), the obvious candidate 
is the hardest jet. Let us call this t Ai ■ That is, 



ta 1 = ^{m 2 i + m\) 



(3) 



where m\ is mass of the hardest jet, (largest energy), and 
m\ is the mass of all the particles not clustered into the 
hardest jet. A quick study of the kinematics shows that 
the t Ai distribution in QCD at order a s is exactly equal 
to the r = 1 — T distribution, with no i?-depcndcncc. 
That is, 
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where the ■ ■ • arc regular in r Al ■ 

Now we turn to the calculation of asymmetric thrust in 
SCET. Note that r A — > forces dijet kinematics indepen- 
dent of which jet is chosen to be the primary one. Thus, 
the effective theory should agree with QCD in the singu- 
lar region for any R. In fact, if we are interested in the 
limit that t a — > 0, then the size of the jet is irrelevant 
to the singularities associated with collinear radiation. 
For this reason, we can use the inclusive jet function. 
Indeed, the use of the inclusive jet function for the hemi- 
sphere mass case has to be justified using the same logic. 
In contrast, for the soft function, the i?-dcpcndence is 
important. The calculation of the soft function is prob- 
lematic for t Ai The problem can be traced to the fact 
that the soft function in SCET has no access to the en- 
ergy of the jets. In the region of phase space in which 
the gluon is soft but not within R of either quark, the 
most energetic jet is not well defined. In QCD, there 
is no such ambiguity. We conclude that SCET cannot 
calculate t Ai . 

Instead, one might consider calculating r Aq , where the 
primary jet is taken to be the quark. When the radiation 
is clustered with the quark, we find that the contribution 
to the soft function from radiation inside a cone of size 
R is 
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where the ^-distribution notation can be found in [6[. 
When the radiation is outside the quark's cone, the con- 
tribution from the complimentary region (everything but 
the cone) is simply S^_ R (k, /u). When we convolute these 
two contributions to the soft function with the inclusive 
(R- independent) jet functions, as in Eq. (TTJ), we produce 
the ta distribution in SCET. The result is 
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Comparing to Eq. (Q| we see that the asymmetric thrust 
distribution in SCET has R dependence, while the QCD 
result for t Ai does not, so the distributions do not agree. 
Of course, it is not unreasonable that they do not agree, 
since they are different observables. 

If we repeat the QCD calculation, but always defin- 
ing the primary jet as the quark jet, we find a result 
which agrees exactly with the SCET prediction, Eq. © , 
up to non-singular terms. Unfortunately, this quark-jet 
asymmetric thrust is not infrared safe, and therefore can- 
not be used to measure jet mass. Indeed, testing t a at 
next-to-lcading order by taking the primary jet to be the 
hardest one containing a quark parton, we find negative 
cross sections. 

Instead of attempting to define a quark jet in QCD 
(see [2(| for an example of how this might be done), we 
can get agreement by averaging the jets. Note first that 
the hardest jet asymmetric thrust, which we calculated 
in QCD, is not the average of the quark and anti-quark 
asymmetric thrust in SCET. However, if we calculate 
asymmetric thrust in QCD with the second hardest jet 
as primary, we find 
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Thus the average of the harder and softer jet asymmetric 
thrust distributions in QCD is the same as the average of 
the quark and anti-quark distributions in SCET at order 
a s . While this may not be surprising, it illustrates the 
importance of being careful not to discuss quark jet and 
gluon-jet masses in QCD without a proper infrared safe 
observable definition [2l|, [22[ . 

Now that the observables agree at order a s we can 
check resummation by expanding the rcsummed results 
in SCET to order a 2 s and comparing against the output of 
the numerical QCD calculation using the program event 
2 [U . We look at the observable 
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FIG. 1: Difference between QCD and SCET for the sum FIG - 2: Difference between QCD and SCET for r Al ^f- at 

ta-i af^j- + ta 2 at order a 2 . That the difference vanishes order a 2 for a variety of jet sizes. 

at very small ta (left side), shows that SCET reproduces all 
of the large logs of this distribution up to NNLL. 

jet radiation, up to a contribution in which the jet mass 
is exactly zero). We will consider several different dijct 
mass observables and explore to what extent SCET can 
reproduce the QCD results and perform resummation of 
the large logarithms. 

The first observable we will consider is the doubly dif- 
ferential cross section in the masses mj and m-i of the two 
most energetic jets. As we saw with asymmetric thrust, 
one cannot order the jets by energy in SCET since, af- 
ter matching, the operator only has access to the jets' 
energies through the labels, each of which has the same 
magnitude. Instead in SCET one can calculate the quark 
and anti-quark jet mass, but these are not infrared safe 
in QCD. As with asymmetric thrust, we will see that a 
meaningful comparison can be made between an average 
over the two jets or more simply the sum of the jet masses 
squared. 

First consider the QCD calculation of d cr/dmfdmf- 
The tree level process e + e~ — > qq and its one-loop virtual 
corrections contribute to the 5{m\)8(m\) since the final 
state partons arc masslcss. This contribution is given by 



and compare to the result of adding the quark and anti- 
quark jet asymmetric thrust distributions in SCET. 

Figure [1] shows the difference between the QCD and 
SCET predictions for ta as a function of lnr^ in the 
singular region, separated by color structure. Figure [5] 
shows the distribution for a number of different values of 
R. In all cases, there is excellent convergence at ta — > 0, 
showing that SCET accounts for all the singularities at 
NLO. This demonstrates that SCET can resum all of the 
large logarithms of ta for any R. At very small R, these 
may not be the dominant part of the distribution, due 
to large lni? terms, but we have not attempted to resum 
logs of R for ta ■ 

As observed in fl~7j , ta may be a promising to measure 
the jet mass in W/Z/j+jet events at hadron colliders. It 
can be calculated at the machine threshold and should 
undergo dynamical threshold enhancement making it rel- 
evant at the partonic threshold as well. An alternative to 
invoking dynamical threshold enhancement would be to 
have a jet mass observable which is not sensitive to the 
whole event. For example, one could look at the mass 
of a jet in exclusively 2-jet events. To do this requires 
a veto on additional radiation. So we now turn to the 
calculation of jet mass with a jet veto in SCET. 
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EXCLUSIVE JET MASS: 
VETO 



DUETS WITH A 



In this section, we will consider the exclusive case, 
where jets are clustered using Cambridge/ Aachen, as be- 
fore, but only the two hardest jets are considered. For 
asymmetric thrust, there was one primary jet, with the 
secondary jet containing all the other radiation in the 
event. For an exclusive dijet sample, we have to veto ad- 
ditional jets. Wc do this by imposing an energy cutoff uj 
on the 3rd most energetic jet: E$ < uj (at tree level and 
0(a s ) this is equivalent to vetoing on the total out-of- 



For the real emission process, e + e~ — > qqg, there are non 
trivial contributions to the jet mass distribution. 

The differential cross section for real emission of a 
gluon in d — 4 — 2e dimensions, leading order in a s , 
is given by 
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S=(Pq+Pq) t=(Pq+Pg) U=(Pq+Pg)- (10) 
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FIG. 3: The allowed kinematic region for e + e~ — > qqg in vari- 
ables u,t. In the yellow regions, Yij, the Cambridge/ Aachem 
requirements have been met and partons i,j have been clus- 
tered. The green region corresponds to E3 < cu. In this plot, 
R = 0.4 and lu = 0.15Q. 



The kinematically allowed regions are shown in Fig. (|3|). 
The yellow regions shown in the figure correspond the 
regions in phase space where the C/A clustering re- 
quirements are met. These regions are labeled for 
it 3 — where the subscript corresponds to to the 

pair of partons that have been clustered. In these regions, 
the clustered jet has the most energy and so we assign 
Pi = pi +Pj, whereas P2 is assigned to the momentum of 
the remaining parton. The shaded green region accounts 
for the case when the partons are not clustered so there 
are three massless jets; however, these are still counted 
as dijct events since E3 < to and thus contribute to the 
5{ra 2 )5(m 2 ,) part of the distribution. After performing 
the integrations, the real emission contribution is given 
by 
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The result of adding the tree, virtual, and real emission 
contributions is 
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As we saw in the last section, SCET docs not have 
access to energy ordering of the jets, and the closest ob- 
servable we can calculate is d 2 a/dm q dm^, which we turn 
to now. The derivation of the factorization theorem for 
d 2 a/dm 2 , dm? is almost identical to the derivation of the 
factorization theorem for thrust. For very small m q _ qi 
the jet size R affects the collinear sector only at sub- 
leading power of m q _q. Therefore, as in the hemisphere 
mass case, we can use the inclusive (i?-indcpcndcnt) jet 
function. The factorization formula then reads 
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This soft function is defined as (see ||) 

S R (k q ,k 9 ,cj,^) = ^<0|Y fi V n |X s )<X s |Y n ty fl | ) 

x 6{k q ~n-p q x )S(k q -n-p 9 x )6{u - k° ut ) , (14) 

where n M and n M are lightlike 4- vectors in the q and q jet 
directions, respectively, Y n are soft Wilson lines, X s is the 
soft radiation, p x {Px) the component of its momenta in 
the quark (anti-quark) jet, and fcg ut is the energy the 
hardest jet of size R in the inter-jet region. 

At order a s , the soft function can be written as a prod- 
uct of two in-cone soft functions and an out-of-cone soft 
function: S R = S%(k q , n)S%(kq, v)S° R ut {u;, fj,). The in- 
cone soft function is the same as for asymmetric thrust 
and given in Eq. ([5]). The out-of-cone soft function is 
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with fo(R) given in Eq. (fTTj) . These results agree with 
limits of expressions given in [l3j . 

Combining this soft function with the hard and inclu- 
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sive jet functions, we get 



The terms singular in t u in full QCD at order a s , are 
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We can now compare SCET directly to the doubly dif- 
ferential distribution of the harder and softer jet mass, 
which is the natural quantity to measure experimen- 
tally, and we have calculated in QCD. The coefficient of 
6(m 2 )5(m^) in the SCET distribution matches the coef- 
ficient of 5(m\)8{rn 2 ) from Eq. (fT2"|) with fo(R) appearing 
instead of f u (R)] however, the remaining distribution is 
symmetric in m q <H- mq in SCET, but not symmetric in 
mi <R- m,2 for QCD. So the doubly differential distribu- 
tions arc different. We anticipated this, but it is still 
instructive to see that the mass of the hardest jet is not 
simply related to any projection of the doubly differential 
distribution in SCET. 

Although the hardest jet mass is not reproduced in 
SCET, one can compute the average of the distributions 
of the two jet masses. If one integrates one jet mass 
up to the full kinematic limit, in general there will be 
non-global logarithms. The result is like comparing the 
left-hemisphere mass in QCD and in SCET, which do 
not agree due to non-global logarithms due to extra hard 
emissions within one jet. However, at small R, there is a 
natural cutoff m < Q 2 R. Then we can compute 
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This can be compared to the analogous expression in 
SCET, using m q and m q instead of m\ and mi. The 
two distributions agree, as they should. 

Rather than look at the average jet mass distribution, 
we can instead look at a thrust-like variable, the sum of 
the squares of the jet masses. We call this jet thrust, 
and denote it by t u . By definition, r w = (m 2 + m\)/Q 2 
where events with £3 < uj have been vetoed. At lead- 
ing order, dcr/dr^ can be calculated by integrating the 
doubly differential distribution computed above. 
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When t w <C 1, each jet mass is forced to be small and 
so SCET is valid in this regime. In the effective theory, 
we take t u = {m 2 + m^)/Q 2 and its distribution can be 
obtained via 
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The terms singular in r w in SCET at order a s , are 
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This result is the same as the QCD result with fo(R) 
appearing instead of f u {R)- Thus at LO, SCET repro- 
duces the parts of QCD which arc singular in t u , up to 
power corrections in uj/Q, but including all of the R de- 
pendence. 



4. PREDICTING HIGHER ORDER TERMS 

In the previous sections we explored doubly and singly 
differential jet mass distributions in SCET and in QCD. 
For some observables, we saw that SCET reproduces the 
singular behavior of QCD at order a s . In particular, 
for jet thrust, r w , which we defined as the sum of the 
jet masses squared normalized to the machine energy Q 
with a energy veto E < u on the third hardest jet, we 
found agreement with QCD for any R at order a s , up to 
power corrections in lo/Q. Next, we would like to know 
to what extent this agreement persists to higher order. 
The question can be phrased succinctly as asking which 
of the coefficients Cij in the expansion 



Tui^— = Y] C ij ( R > W ) Q ' S lnJ T ^ 



(22) 



can be predicted in SCET. 

There are two parts to this question. First, if one were 
able to do calculations to all order in SCET, would the 
coefficients come out correct? As long as the entire sin- 
gular region is described by soft and collincar degrees of 
freedom, the answer should be yes. This is true for t w as 
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/ C F 2 (withoutNAE) 

<r CprifTp 








R=0.1,u= 1(T 4 Q 





-10 -8 -6 -4 -2 



FIG. 4: Difference between QCD and SCET for the (a s /4?r) 2 
coefficient of T u £f- for the different color structures. After 
using non-Abelian exponentiation (NAE) the Cf curve should 
go to zero at small r„, up to power corrections in oj. 

long as ui <C Q for any R. So SCET should get all the 
Cij(R,ui) correct up to power corrections in ui/Q. This 
can be seen already at order a s in Eq. pip where the 
SCET distribution agrees with the QCD distribution in 
Eq. (fH?|) exactly in R when terms of order uj/Q or smaller 
are dropped. 

The second part is whether the coefficients dj(R,uj) 
with i > 1 can be predicted using factorization. That 
is, which of these coefficients can we get without ac- 
tually doing the explicit loop integrals? For example, 
non-global logarithms appear in the doubly differential 
hemisphere mass distribution. These are reproduced in 
SCET [13, [ll| , but it is so far unclear how to reproduce 
them without doing calculations to all orders in pertur- 
bation theory. We would like to know what can and 
cannot be predicted, in what limits. We would also like 
to know to what extent terms which we cannot predict 
or calculate are actually important. 

Since the anomalous dimensions of the jet and hard 
functions are known to 3-loops, the /i-indcpcndcncc of 
the t u distribution determines the anomalous dimension 
of the soft function to 3-loops as well. Note that since 
the jet and hard functions are completely independent 
of R and w, none of the anomalous dimensions have any 
dependence on R or ui at all. With the finite part of 
the 1-loop soft function known, this should naively al- 
low us to resum the next-to-next-to- leading logs (NNLL). 
That would give us all the terms down to 1 /t u in the a 2 
distribution. This is certainly the case for hemispheres 
(R = |). However, for smaller R these anomalous di- 
mensions are not enough to predict the NNLL, or even 
NLL distribution. 

The failure of SCET to describe at NLL can be 
seen already in Eq. (j2"Tj) : at NLL, the 1-loop anomalous 
dimensions should determine the a s /r w terms. But since 
the jet and hard anomalous dimensions lack R depen- 
dence, there is no way they can predict the a s In jz^/t^ 
piece. Thus, without insights into the soft function be- 
yond the original factorization, SCET cannot even resum 
t„ to NLL. 



To proceed, we first of all observe that the troublesome 
term has a coefficient Cf- Due to the non-Abelian ex- 
ponentiation (NAE) theorems [H, [H| the a n 8 C F terms 
in the soft function are determined exactly by the 1- 
loop result. Thus we can resum at least this much in- 
dependence by simply exponentiating the NLO soft func- 
tion (in position space). Doing so predicts an NLL term 
Cpa 2 In 2 t/t, as well as a 2 In t^/t^ and a 2 /r w NNLL 
terms. 

To see how well SCET is doing, we will explore the or- 
der a 2 s distribution of t u . This can be calculated in QCD 
numerically using the program event 2 [25[. In Figure 0] 
we show difference between QCD and SCET predictions 
at order a 2 for the different color structures. The Cp 
curve before and after using non-Abelian exponentiation 
are shown. We see that after NAE, the Cp curve goes 
to zero as t u — > 0, while the CfCa and CprifTp color 
structures are clearly missing \t\t u Q/uj pieces, which af- 
fects the slope of these plots. Note that the curves cross 
at roughly t u ~ cu = 10~ 4 ~ e~ 9 , indicating a missing 
In ui /t u piece. The approach of the Cf term to zero gives 
us a benchmark for where we should expect the curves 
to vanish, and how large the uj/Q power corrections are. 

A. Refactorization 

To further progress, we will now argue that when 

u<k Ll k R and R < 1 (23) 
the soft function factorizes into four parts 

®S f (k L ,k R ,uj), (24) 

Here, £g> indicates a convolution. Similar refactorization 
formulas have appeared before [ill 1111 (and after [26j). 
but of a more general nature. Without specifying the 
regime in which the refactorization is valid, and more 
importantly, properties of the objects appearing in the 
refactorization, there is no content in Eq. (j2~4"l) . 

To see that Eq. by itself is contentless, recall that 
the hemisphere soft function can be written as 

<Shemi(&L,fcR,M) = n(fc if fJ,)U(k R , /J,) !g) Sf(k L , k R ) , 

(25) 

where n(fc, /i) are the soft RG factors which compensate 
for the RG evolution of the jet and hard functions. This 
form is forced by the hemisphere factorization formula in 
SCET [Ell, But S hemi = 5^1, so we can find 

a trivial solution for f24|) by setting S l R (k,fi) = U(k, (j,) 
and ^ ut = 1. 

We propose that S R and S R ut can be thought of as soft 
functions in their own right, with their own anomalous 
dimensions. Only has a log /i in its anomalous di- 
mension, since the soft-collincar singularities are within 
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the jets. We conjecture that the regular anomalous di- 
mensions split as 



7s Ut = -r cusp ln 



R 



7A 



1 - R 
R 

7s = r cus P In Y~7~r + 7s - 7A 



(26) 



where T cusp is the cusp anomalous dimension. In particu- 
lar, all of the i?-dependence in the anomalous dimensions 
in proportional to r cusp . 

We will argue that 7a = as well, which is consis- 
tent with the R = \ limit, although since our refactor- 
ization argument is valid at small R, this is not really 
justified. We also set Sf(kL, kn, uj) = 1, which is to say 
we ignore non-global logs. We will see that the refac- 
torization is still in good agreement with full QCD in- 
dependent of the non-global log issue. With this conjec- 
ture, the refactorization allows us expand around small 
R to predict non-trivial parts of the t u distribution, for 
u/Q < ^ R <^ 1, which we will compare our predic- 
tion to the fixed order result. 

A complete proof of refactorization (or more precisely, 
why 5^ ut (a;,/i) ^ 1) would ideally involve operator defi- 
nitions of the components and S^ ut . That way we can 
ascribe anomalous dimensions to the components, which 
can, at least in principle, be calculated order-by-order in 
perturbation theory. However, it may not be possible, 
or even necessary, to have operator definitions for refac- 
torization. For example, in the hemisphere case, we do 
not have operator definitions of the components of Shemj. 
Instead, Tl(k, /i) is defined to compensate the hard and 
jet function evolution, and Sf is defined as whatever is 
left over. We aim at this point for the more modest goal 
of having a similar algorithmic definition of the compo- 
nents of Sr, leaving a more complete non-pcrturbative 
understanding for future work. 

In Section 14 Bl we give a heuristic argument for why 
the soft function should refactorize at small R (and 
w ^ We present phenomenological evidence in 

support of this refactorization conjecture in Section |4 CI 



B. Heuristic argument for refactorization 

To justify the refactorization, we will show that the 
kinematic restriction small R and uj < k^^kn allows us 
to apply soft-collinear factorization [29(. The argument 
we present here is similar to one used in [30| for when 
two subjects become collinear inside a larger jet. 

To begin the argument, let us define the "left" jet to be 
in the n M direction and the "right" jet to be in the di- 
rection. A particle within the left jet is kinematically re- 
stricted to have fcz, = fc + < j^^k~, where k~ = n-k and 
k + = n ■ k. With the on-shell condition, the transverse 



momentum then scales like k i 



Vk+k^ ~ \lt=Ti k ~ 



So, when R <C 1, the jet is forced to have collinear scal- 
ing (k~, k + , fcj_) ~ %■(!, R, y/R). Radiation in the right 



jet is restricted analogously. Outside the jet, the soft mo- 
mentum is only required to have fc out ~ fci(l, 1, 1). Thus, 
there are three sectors, |A" ou t), which has ultra-soft scal- 
ing and \X£) and \XP) which are soft, with respect to 
the original collinear momentum, but collinear with re- 
spect to the soft radiation outside the jets. We say these 
have soft-collinear scaling. The relevance of soft-collinear 
modes was also noticed in [30| in a different context. 

It is interesting to observe the soft-collinear modes 
are formally harder than the original soft modes: in 
terms of r u ~ k^/Q, the soft collinear modes scale like 
Qt u (^, 1, ^t=) while the original soft modes scaled like 
Qt u (1,1, 1). However, the soft-collinear modes are still 
parametrically softer than the original collinear modes, 
which scale like Qt u (-£-,1,-^}=), since t u <C R- Thus, 
both the soft-collinear and ultrasoft modes have eikonal 
interactions with the original Wilson lines, and hence are 
present in the original soft function. 

Now consider the original Wilson line 



Y n (x) = Pexp 



i ft 



ds n ■ A(x + sn) 



(27) 



The gauge fields in this Wilson line includes anything 
with Eikonal interactions with the original collinear 
fields, which includes both the ultra-soft and soft- 
collinear modes. Thus, we can write = A™ 8 + A™ 5 + 
A" s . The soft-collinear modes and ultra-soft modes intcr- 
act with each other the same way that ordinary collinear 
and soft modes interact. In any situation in which modes 
have collinear and ultra-soft scaling, the interactions be- 
tween the soft and collinear modes can be removed from 
the Lagrangian through a BPS field redefinition [2j| : 

A ns -+ (Y^A ns Y™, A ns -+ (Y^A^Y™. (28) 

where Y" s is the same as Y n but with only the ultra-soft 
gauge fields involved. After this redefinition, there are 
no-longer interactions between soft-collinear and ultra- 
soft modes. Thus, the original Wilson line operator now 
separates into 



yty 

1 n 1 n 



(Y£ c )\Yry(YD(Y° c ) 



(29) 



where Y£ c is a soft-collinear Wilson line, again given by 
Eq. ((27)) . but now in terms of decoupled soft-collinear 
fields. This separation of the soft Wilson line into soft 
and soft-collinear Wilson lines suggests the refactorized 
form in Eq. ([2"g]l. 

Note that this derivation is not a complete proof of the 
refactorization. We have not given operator definitions 
for and iS^ ut , although clearly they should be related 
to matrix elements of the soft-collinear and ultrasoft Wil- 
son lines respectively. This section should be viewed as 
an argument for why the soft function might simplify as 
R — > 0. Such simplifications would show up in the calcu- 
lations of the exact higher-order soft function, which is 
not yet known. If the reader does not find this argument 
convincing, he should ignore this section and consider the 
phenomenological evidence for refactorization at small R 
which we give below. 
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C. Predictions from refactorization 

To solve the RGE for the soft function, we first ex- 
tract the anomalous dimensions of the two parts from 
their NLO expressions. The anomalous dimension of 
5^ ut (w, fi) has no \i dependence, consistent with there be- 
ing no soft-collinear radiation in the middle region. The 
/i dependence of the anomalous dimension of S^(k, is 
the same as for the thrust soft function. We can also read 
off the regular anomalous dimensions at NLO 

IT —Itlpln^ (30) 

7^ = 75 + rW p i nT A_ 

where Tculp = CfC( 3 /tt is the one- loop cusp anomalous 
dimension. Generalizing this to higher orders leads to 
our conjecture Eq.fJSB]) with 7a = 0. As discussed ear- 
lier, we are not interested in non-global logs, so we set 
Sf{kL,kn 1 oj) = 1 for simplicity. 
In summary, our conjecture is 

S a (k Lt k a ,LJ,n) = Sj?(fc L ,M)^(W)Sr(^M) (31) 

with Sft(k,fi) having the \n(fi/k) part of the anomalous 
dimension, and the regular anomalous dimensions split- 
ting as 7g ut = -r cusp In and 7^ = 7s + r cusp In 

Although we have not proven the soft anomalous di- 
mension splits up in this way at higher order, we do know 
that the R dependence must exactly cancel to all orders, 
which implies it should be universal, naturally suggest- 
ing the cusp anomalous dimension. As a check, we also 
know that as R — > i, 7i n — > 75 and 7 ou t — > 0, since this 
is the hemisphere case. At order a 2 s: this r cusp ansatz 
contributes terms like Fi In lnr w to the singular t u 
distribution, which we confirm exist by numerical com- 
parison to QCD (see below). We do not have a proof 
for the r cusp ansatz or a strong argument about why the 
regular anomalous dimension cannot split up differently 
at two loops and beyond. 

We also note that the connection between the lni? 
terms and the cusp anomalous dimension was pointed 
out in [l3|, [l5|. In fact, part of the reason we expect that 



our r cusp ansatz could hold is related to observations 
about the structure of the divergences in the relevant 
Feynman diagrams which generate the 1-loop anomalous 
dimensions, first calculated in these papers. However, 
without the refactorization at small R, there is no pre- 
dictive power in this observation. Indeed, generically, 
there can be an additional anomalous dimension 7a, as 
indicated in Eq. (|26p . and 7a can also depend on R. 
We are arguing that all the singular dependence on R 
in the anomalous dimension should be given by the cusp 
anomalous dimension, so that at small R, this 7a is it- 
independent (and vanishes at one-loop). We cannot rule 
out such a function which vanishes at R = i and goes to 
a constant starting at two loops. The small R regime in 
not mentioned at all in [ill, El (although it is discussed 
in a paper (3~lj on non-global logarithms). 

As a consequence of refactorization, we can choose sep- 
arate soft scales for the in-jet and out-of-jet regions, /x s 
and Evolving to [i u is straightforward, since 
there are no scales in its anomalous dimension. Thus 

S °* { "^ )= \J=r) ST^,^)- (32) 

The function Ay (y, (J.) is defined in [l8| . It follows that 
choosing fi u = u allows us to resum logs of u), since there 
is no other scale in this soft function. 

The soft function is similar to the thrust soft func- 
tion up to i?-dependcncc. Its evolution equation is simi- 
larly solved in Laplace space [IH . The projection relevant 
for t w adds the /cl and kn momentum is just the thrust 
soft function times a new independent factor. 

While there may be non-global logs contained in the 
function Sf(kL,kn,ui) in Eq. (f24f. for simplicity, we have 
set Sf = 1. Non-global logs are not the subject of this 
paper. To the extent that our refactorization agrees with 
the NLO, it is partly because non- global logs arc a nu- 
merically small effect in the regime we consider. We refer 
the reader to [26| , which appeared after the first version 
of our paper appeared, and addresses non-global logs for 
jet thrust and related observables. 

Combining the full soft function with the hard and 
inclusive jet functions gives the resummcd t u distribution 
in SCET 



1 

(T{) 



'da 

dTui 



= exp 



SCET 



4S(fj,h,fJ>j)+4:S(fj >s ,fj,j)- 2A H (ph, fJ-s)+^Aj(fij,fi s ) 



R 



1 - R 



91 
/'/, 



-2A T {n h ,tij) 



xff(Q 2 ,^)^ ut (^ 



J[ In— 2- 



2£ (d^/z, 



1 ( tM 

Us 



71 e ~lEV 



(33) 



Note that this expression is explicitly independent of scales fixes the unknown parts of the fixed order soft 
(i. Demanding independence of the various matching functions. Both are of the general form of fixed order 
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FIG. 5: The difference between coefficient of in dcr/dhnvj 
in full QCD and in SCET for R = 0.1 and u = 0.0001Q after 
refactorization, but not including the Fi piece. 



FIG. 7: The difference between coefficient of a 2 a in dcr/dhnvj 
in full QCD and in SCET for R = 0.3 and oj = 0.01Q after 
refactorization and including the Fi piece. 
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FIG. 6: The difference between coefficient of o? a in dcr/dlmvj 
in full QCD and in SCET for R = 0.1 and u = 0.0001Q after 
refactorization and including the Ti piece. 



expansions in RG-improved perturbation theory, given 
by h(L,n) in Eq. (55) of [H]. In this case S% ut (w , fx u ) = 
/i(21n^,/i s ) with all the r cusp terms set to zero, since 
there are no Sudakov logs outside of the jets, and sub 
stituting 7# - 



7 out 



rcusphij^j and 



-7g ut with 

2C F In 2 j=ji + fo(R). The in-cone soft func 



cf -> c° u 

tion has double logs, and is given by (L, /i) = /i(2L, //) 
with 7 ff -> -7^, 7^ = 7h - 27j + r^spln^, and 

cf ci 11 = -Cf(7t 2 + 2 In 2 y^j). Since all the anoma- 
lous dimensions are known up to 3-loops, the first un- 
known terms are the 2-loop constants in the soft func- 
tions. Thus, the refactorization formula and the r cusp 
ansatz, that the 2-loop out-of-cone soft function anoma- 
lous dimension is proportional to Ti, gives us enough for 
NNLL resummation (again, up to NGLs). 

To check our result, we can expand our result to order 
a 2 and compare the resulting distributions to full QCD 
using event 2. The results should agree up to power 
corrections in t u , lo/Q and R and up to NGLs. Figure [5] 
shows the difference between SCET and QCD as a func- 
tion of lnr w for R = 0.1 and uj/Q = 0.0001 using refac- 



torization (but setting the 2-loop 5 out anomalous dimen- 
sion to 0). We see that now both the CaCf and C^n/Tp 
color structures go to zero slope, similar to what was seen 
already in Figure |4] for the C\ color structure (recall that 
the Cp part of the 2-loop soft function is known exactly 
from non-Abelian exponentiation). This indicates that 
we are getting the In t u /t u terms mostly correct. We 
then show in Figure [6] the further improvement gained 
by including the contribution from ri in the 2-loop soft 
function anomalous dimension. That the curves now go 
to zero is a highly non-trivial check. This indicates that 
we are getting also the l/r w terms mostly correct. The 
refactorization is easiest to confirm in the region with 
t u ~ to, where the NGLs are necessarily small. 

As discussed above, the refactorization we presented is 
valid only for small cj < t^Q and small r w <C R <C 1, 
but we expect it to be relevant even for moderate R. 
In Figure [7| we show the case with R = 0.3. One can 
see clearly that the refactorization is phenomenologically 
important even for R which is not terribly small. 



5. CONCLUSIONS 

In this paper we have considered a number of impor- 
tant issues which will be relevant to producing accurate 
distributions of jet masses and other jet substructure ob- 
servables at colliders. We first considered the possibil- 
ity of an inclusive jet mass measurement, such as the 
asymmetric thrust observable. In this case, we were able 
to resum logarithms of the jet mass to next-to-next-to- 
lcading logarithmic accuracy. These distributions, when 
expanded to order a 2 showed excellent agreement with 
numerical calculations in full QCD in the threshold re- 
gion. 

At hadron colliders, it will be difficult to use an in- 
clusive jet mass observable. Instead, one must restrict 
consideration to the hardest jet or the hardest pair of 
jets in the event. In this case, it is important to veto 
additional jets, which introduces a veto scale u). Then 
there are at least three parameters, oj/Q, where Q is 
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some hard energy scale, the jet shape r and the jet size 
R. The goal in resumming the jet shape would be to 
predict coefficients of a\ In 7 r without doing the full cal- 
culation to ith order in perturbation theory. In general, 
this is difficult, but we argued that there is a refactoriza- 
tion which simplifies the calculation for the observable 
jet thrust, t u . Jet thrust is defined as the sum of the 
jet masses squared normalized to the center of mass en- 
ergy Q, after a vetoing the energy uj of the third hardest 
jet. We argue that the refactorization holds in the limit 
w/Q < < < 1. We presented kinematic scaling 
arguments, using the language of Soft-Collinear Effective 
Theory, to justify the refactorization. 

In order to test the refactorization, we compared the 
predictions to the exact NLO event shape distributions in 
full QCD. This provided very strong numerical evidence 
that the refactorization correctly predicts the dominant 
part of the jet mass distributions. We further argued that 
the i?-dependence may be universal, and therefore may 
depend on the cusp anomalous dimension. This T cusp 
ansatz was also confirmed by comparison to full QCD. 
The predictions from refactorization are independent of 
the separate, and important, issues of non-global logs. 
NGLs arc apparently a numerically small correction to 
our results, but arc not the focus of this paper. Extrap- 
olating away from R = our result still provides good 
agreement with full QCD. Thus the small R expansion 
seems like a productive direction for further investiga- 
tion. 



The improved agreement with full QCD after refac- 
torization is impressive, but still not completely well un- 
derstood. It would be helpful to have the full soft func- 
tion calculated, along the lines of [l(| [H| ■> to understand 
exactly which terms are reproduced and which are not. 
It would also be helpful to have gauge-invariant opera- 
tor definitions of the functions into which the soft func- 
tion factorizes, so that they can be computed directly. 
Progress along these lines has already appeared, in [30{ . 
Most importantly, it will be critical to construct jet mass 
observables for hadron colliders which can be computed 
and resummed to next-to-next-to-leading logarithmic ac- 
curacy so that direct comparisons to collider data can be 
made. 
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Appendix A: NNLO expansion of SCET results 



In this appendix, we give the t w distribution in SCET to order a 2 s . 
Eq. (|2~Tj) . At order a 2 , we can write 



The SCET distribution at order a s is given in 



da 

dTu, 



SCET 



-( 



' NAE 



CfCa 



cCa _r_ f<^_ 

Jo + iri 



C A 



f0 f +/r/] +C F /? /Rcfact} (Al) 



where /?o = ^-C-a ~ ^ n fTF is the leading /3-function coefficient in QCD. 

The parts coming from the 1-loop result and the original SCET factorization formula for r w are 



f° F =321n 3 r w + ln 2 r w 



72 - 96 In ■ 



R 



1-R 



lnr„ 



32Li 2 (l - R) - 32Li 2 i? + In | 64 In ^ 4 



24Li 2 (l - R) - 24Li 2 (i?) + In — — ( 48 In ^ + 24 lni? + 24tt 2 - 16 ) + I6C3 - 8tt 2 + 9 (A2) 

1 — R \ Q 



f° A =441n z r w 



ln- 



352 



In. 



R 
l~R, 



r 



338 



176, 2u, 

In — In — 

3 Q 1 



R 



R 



— Li 9 (l - R) Li 2 i? H In - 

3 V ' 3 3 1 



R 



R 



lni?+24C 3 -57 (A3) 
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f lf = -161n 2 r w + ln< 



128 , R 
In ■ 



3 1 - R 9 



The part coming from non-Abclian exponentiation is 

= 64 In 2 — — lnr w +321n— — ( 2Li 2 i? - 2 In — — In — - In — — ln(4i?) + lnfl - R) lni? - — ^ (A5) 
JNAE 1-R 1-R\ l-R Q l-R K 1 y 1 12 J y ' 

The part from the soft-refactorization is 

/Refact = 16 In ^ In — + ALh(R) + ln 2 (l — R) + 2 ln(l - R) In R - In 2 R - 1_. (A6) 
2a; 1 — it 3 

The parts from the I\ ansatz are 

/ r ^ = 5(67 - 37r 2 )ln-^-, tf^-i^ln-J* (A7) 

Adding these pieces, the complete CfCa color structure at order a 2 is relatively simple 

, Ci ... a /8tt 2 176, i? 338\, 536, i? 44 , 2 R 8tt 2 , i? n _ rw 
jOa = 441n 2 + m m H j 1 j 2. j 24 > _ 57 

J V 3 3 1 - # 9/ 9 1 - i? 3 1 - i? 3 1 - ii s 

(A8) 

Note the w dependence has completely dropped out of the sum for the CfCa (and the Cf^/Tf) color structure order 
a 2 . This follows simply from the refactorization. There is u dependence in the C F terms. 

By explicit calculation, we can also work out the non-global log in the CfCa color structure. Rather than use the 
Cambridge/ Aachen algorithm, for simplicity, we calculate the contribution by just including all radiation within R of 
the jet direction. This produces a non-global log of the form 

c c A ( 16tt 2 , onT: R 2 \, (t u Q 



&& = i-— +32Li2 (T3^J 1 Hi?)- (A9) 



This vanishes at R — | and goes to — In at R = 
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